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ABSTRACT 

We present a novel approach to quantifying the morphology of Cosmic Microwave 
Background (CMB) anisotropy maps. As morphological descriptors, we use shape 
parameters known as Minkowski functionals. Using the mathematical framework pro- 
vided by the theory of integral geometry on arbitrary curved supports, we point out 
the differences to their characterization and interpretation in the case of flat space. 
With restrictions of real data - such as pixelization and incomplete sky coverage, to 
mention just a few - in mind, we derive and test unbiased estimators for all Minkowski 
functionals. Various examples, among them the analysis of the four-year COBE DMR 
data, illustrate the application of our method. 

Key words: methods: numerical; methods: statistical; cosmic microwave background 



INTRODUCTION 



The oldest signal accessible to mank ind is the Cosmic 
Microw ave Background discovered by Penzias & Wilson 



ical defects such as strings or te xtures arising from phase 
transitions in the early Universe (stebbins 1988). 



(1965). Consisting of photons that have been free-streaming 



since the Universe was only 300,000 years old, the Cosmic 
Microwave Background (CMB) provides valuable informa- 
tion on the early history of our Universe. Above all its 
anisotropies mirror the matter fluctuations at the epoch of 
recombination (at redshift z « 1, 100) and with them the 
seeds of the Wgp-gralp structure today 



In order to measure the morphology of Cosmic Mi- 
crowave Background anisotropies, the Euler characteristic, 
or e quivalently the genus, was suggeste d as long as a decade 
ago (Coles & Barrow 1987, Coles 1988). Even to date, most 
applications are confined to genus statistics, alt hough an 
early theoretical study by Gott III et al. (1990) also con- 



siders the boundary length, but failed to come up with a 
subsequent analysis of data. The analysis of the fir st-year 
COBE DMR data using genus statistics was done by 



Smoot 



et al. (1994); their work also contains a thorough discussion 



Various methods of statistical analysis have been used 
on Cosmic Microwave Background maps. Amo ng them arc 
the two- and three-point cor relation function ( Hinshaw et 
al. 1996 L Hinshaw et al. 1995), the power spectrum flGorski 
rt"a!"T996|) , skewness and kurtosis (Luo fc Schramm 1993|), 
multifr; £tals_( Pompilio et al. 1995 ), and the cxtrema corre- 
lation function (Kogut et al. 1993). Another promising ap- 



of the performance of the method compared to other mea- 
sures of non-Gaussianity. Further applicat ions of topolog ical 
methods on CMB anisotropies come from Torres (1994) and 



Torres et al. (1995), and gen us calculations of th e fo ur-year 



COBE DMR data are due to |Colley et al. (1996)| and |Kogut 



et al. (1996) 



proach is the investigation of the morphology of hot and cold 
spots. Being complementary to the traditional approach via 
the hierarchy of correlation functions, it provides al ternative 
method s for determining cosmological parameters ( Torres et 
al. 199E ). But above all, morphological statistics incorporate 



The genus can be pl aced in the wider framework of the 
Minkowski functionals (Minkowski 1903), by natural and 



correlalion functions of arbitrary order. Hence they are sen- 



compelling mathematical considerations. Originally intro- 
duced to tackle problems of stochastic geometry, this family 
of morphological descriptors sub sequently set o ff the devel- 
opme nt of integral geomet ry (sec Blaschke 1936 or Hadwiger 
1957 for early works, and 



Schneider 1993 



sitive to signatures of non-Gaussianity in the temperature 
fluctuations, which would indicate the presence of topolog- 



for a comprehen- 
sive overview). Recently, the Minkowski functionals have 
been introduced into cosmology as descriptors for the mor- 



© 1997 RAS 



2 Jens Schmalzing and Krzysztof M. Gorski 



d 



1 



a linear combination of only d + 1 functionals; these are the 



Vo length 
Vi X 

v 2 
v 3 



area, 
circumference 
X 



volume 
surface area 
total mean curvature 
X 



Table 1. Some of the d + 1 Minkowski functionals in d— 
dimensional Euclidean space may be interpreted as familiar geo- 
metric quantities (apart from numerical factors). This table sum- 
marizes the geometric interpretations of all Minkowski functionals 
for one, two and three dimensions. The symbol \ denotes the Eu- 
ler characteristic, a purely topological quantity; it measures the 
connectivity of a set, being equal to the number of parts minus 
the number of holes in two dimension. 



phological properties of large-scale structure by Mecke et al | 
(1994)[J While their original approach uses a Boolean grain 



model applicabl e to the analysis of point sets, Schmalzing & 



Bucher ; (1997) consider excursion sets and isodensity con- 
tours of smoothed random fields. Applications have so far 
been restricted to the morphome try of large-scale stru cture 
in redshift catalogue s of galaxies (jKerscher et al. 199"7a ) and 
clusters of galaxies (Kerscher et al. 1997t). 

The promising results in large-scale structure analysis 
motivate the application of Minkowski functionals to Cosmic 
Microwave Background sky maps. Since all-sky maps live 
on a curved support, some formal obstacles will be encoun- 
tered, but the underlying concepts remain the same, an d 
the central formulae are easily generalized (Bantalo 1976). 



In applications to data, care must be taken to remove the 
effects of usually incomplete sky coverage, while retaining as 
much information as possible at the same time. 

Our article is organized as follows. Section ^ summarizes 
the framework of integral geometry first in flat space, and 
then on spaces of non-zero but constant curvature, with spe- 
cial regard to the different interpretations. Also, we express 
all Minkowski functionals of the excursion set of a smooth 
random field as integrals over purely local invariants formed 
from derivatives. In Section ^ we use these integration for- 
mulae to construct estimators for the numerical evaluation 
of the Minkowski functionals for a pixelized CMB sky map. 
The important problem of incomplete sky coverage is ad- 
dressed, and we find prescriptions for boundary corrected, 
unbiased estimators, even for smoothed data. Section ^ is de- 
voted to examples, among them a study of noise reduction 
through Gaussian filtering, and a morphological analysis of 
th e maps constructed fro m the COBE DMR four-year data 
by Bennett et al. (1996)| . Finally, we summarize, draw our 
conclusions and provide an outlook in Section H. Two ap- 
pendices further illuminate the mathematical aspects of this 
paper by giving detailed derivations of important formulae. 



2 THEORY 

2.1 Integral geometry 

Let us first introduce integral geometry in flat space, or, to 
be more precise, in d-dimensional Euclidean space E d . We 
wish to characterize the morphology of a suitable set Q C 



E . Hadwiger's Theorem ( Hadwiger 1957 ) states that under 



to d. If the set Q has a smooth boundary dQ, its Minkowski 
functionals - except for the d-dimensional volume Vo, which 
is of course calculated by volume integra tion - are given by 
simple surface integrals (Schneider 1978). So altogether we 
havd^ 



V (Q) = J dv, 



Wj-l \ ■ J dQ 



ds&j (ki ... Kd-l) ■ 



(2) 



Here dv and ds denote the volume element in E and the 
surface element on Q, respectively, /ti to Kd-i are the bound- 
ary's d— 1 principal curvatures, and &j is the jth elementary 
symmetric function defined by the polynomial expansion 

d— 1 d 

\\i x + K = /. x d ~ j &j («i • - • k<i-x) ; (3) 
»=i j=i 

hence ©i = 1, S2 = «x + • • • + Kd-i, and so on up to &d = 
Ki . . . Kd-i- Table |l| summarizes geometric interpretations of 
the Minkowski functionals in one, two and three dimensions. 

2.2 Spaces of constant curvature 

Let us now consider the d-dimensional space of constant 
curvature kK. The sign k equals +1, or —1, for the spher- 
ical space S d , the Euclidean space E d and the hyperbolic 
space HI d , respectively. K is a positive constant of dimen- 
sion [Length] -2 , hence its inverse squa re root K~ x ^ 2 c an be 
interpreted as the radius of curvature, santalo (1976) shows 



how to obtain an integral geometry on such spaces. Curva- 
ture integrals as in Equation (^) can still be defined, if care 
is taken to use the geodesic curvatures Ki. In the following, 
we will call these quantities the Minkowski functionals in 
curved spaces. 

However, some of the geometric interpretations are al- 
tered with respect to the flat case. While in flat space the 
curvature integral Vd(Q) is equal to the Euler characteris- 
tic x(Q)y curved spaces require a generalized Gauss-Bon net 
Theorem proved for arbitr ary Riemannian m anifolds by Al- 



lendoerfer fc Weil (1943)| and |Chern (1944J The theorem 
states that the Euler characteristic is a linear combination 
of all Minkowski functionals as defined by Equation (H), 



x(Q)=X>VKQ), 
3=0 



with the coefficients Cj given by 

d \ 2(kK) (d - j)/2 J . 

. 1 — : 11 d — J even, 

J ) Ud-j 

if d - j odd. 



(4) 



(•>) 



* We use uij to denote the surface area of the j'-dimensional unit 
sphere. Some special values are u>o — 2, w\ = 2n, u>2 = 47r, while 
in general 



a few simple requirements, any morphological descriptor is 
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Note that from the point of view of Hadwiger's theorem, 
which is also valid on curved spaces, all linear combinations 
of Minkowski functionals are equally suitable as morpholog- 
ical descriptors, so one may both use the integrated curva- 
ture Vd and the Euler characteristic \ as the last Minkowski 
functional. In the following, we will consider both quanti- 
ties, because the integrated geodesic curvature is easier to 
calculate, and the Euler characteristic is easier to interpret. 
Obviously, in the case of Euclidean space E d , k = and all 
coefficients apart from Cd vanish, so \ = Vd and the original 
Gauss-Bonnet theorem is recovered. 



2.3 Two— dimensional unit sphere 

We now focus attention on the supporting space for CMB 
sky maps, the sphere § 2 of radius 7?. The parameters intro- 
duced in the previous section now take the values d = 2 for 
the dimension, k = +1 for the curvature sign, and K — R~ 2 
for the absolute value of the curvature. 

Rewriting the definition in Equation (^), we obtain the 
Minkowski functionals for a set Q C § 2 with smooth bound- 
ary dQ by 



V (Q) = J da, 
Q 

Vi(Q) ~ / 'i'- V 2 (Q) = i / da,-. ((.) 



delta function, and the appropriate Jacobian. 



Vi{Qu) = -g J d£ ' 
V 2 (Q„) = ^ j d£ K - 



da S(u — v) \ grad u\ 



da 8{u — v) | grad u\ —k. 

Z7T 



(10) 



Since the integrands can now be written as second-order 
invariants (see Appendix ^ for a detailed calculation of the 
geodesic curvature n), we have succeeded in expressing all 
Minkowski functionals as surface integrals over the whole 
sphere S 2 , 

Vj{Qu) = j dalj, (11) 

s 2 

with integrands X, depending solely on the threshold v, the 
field value u and its first- and second-order covariant deriva- 
tives. In summary, 

To = B(u - u), 



Ti = -5{u- u)^^ + u%, 

^ 1 .,/ , 2u ; lU;2M;12 — M 2 lM;22 — W 2 2 U ; il 
J 2 = — d(u - V) ■ ■ — ^ — — ' 



(12) 



In the following, we will use the surface densities of the Min- 
kowski functionals, that is divide by the area of S 2 . 



AtyR 2 1 3 



(13) 



HQ 



where da and dl denote the surface element of S 2 and the 
line element along dQ, respectively. Being a linear object, 
the boundary dQ has only one geodesic curvature k. 

Using the generalized Gauss-Bonnet Theorem in Equa- 
tion (0) with the coefficients for two dimensions substituted, 
we can calculate the Euler characteristic x(Q) from the Min- 
kowski functionals via 



x(Q) = V2(Q) + 



2-kR 2 



V (Q). 



(7) 



Note that by inserting the definitions from Equation (|6|), 
this formula reproduces the ordinary Gauss-Bonnet The- 
orem for surfaces with a smooth boundary embedded in 
three-dimensional flat space. 

Let us now consider a smooth scalar field m(x) on S 2 , 
for example the temperature anisotropics of the Microwave 
sky. We wish to calculate the Minkowski functionals of the 
excursion set Q v over a given threshold v, defined by 



?„ = { x e S 2 u(x) > v } 



(8) 



The zeroth Minkowski functional Vo, i.e. the area, can be 
evaluated by integration of a Heaviside step function over 
the whole sphere 



V Q {Q V ) 



da Q(u 



(9) 



s 2 



2.4 Expectation values for a Gaussian random 
field 

Minkowski functionals and other geometric characte ristics 
of Gau ssian random fields are extensively studied by Adlei 



(1981). Analytical expressions for the average Minkowski 



functionals of a Gaus sian random fie ld in arbitrary dimen- 
sions were derived by Tomita (1986); in the special case of 
two dimensions, the results for the isodensity contour at 
threshold v areQ 



/ \ 1 It/' v 

VnM — $ — ==■ 

0K ' 2 2 



-1/2 



exp 



2(7 



(14) 



2tt 3 / 



T V — 11 I (V — ll)' 

exp 1 



2(7 



2a 



Note that these expressions contain only three parame- 
ters, namely fi, a, and r. All three are easily estimated from 
a given realization of the Gaussian random field, by taking 
averages of the field itself, its square, and the sum of its 
squared derivatives; then 



A 4 = <") 



/ 2\ 2 

a = (u ) — fi 



T = - (u-iU-i) 



(15) 



The other Minkowski functionals are actually defined by line 
integrals along the isodensity contour in Equation (^|), but 
they can be transformed to surface integrals by inserting a 
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derived from a po wer-law spectrum P(k) oc k n by Bond & 



0.00-1 



0.008 




Figure 1. The angular power spectrum shown in this plot is 
used for all tests of the estimators presented in this section. It was 
chosen as a rough model for the DMR maps while remaining an- 
alytically tractable. The complete power spectrum, indicated by 
a solid line, is the sum of two contributions, namely a Harrison— 
Zeldovich spectrum normalized to Cio = (7/^K) 2 , convolved with 
a 7° FWHM beam to mimick cosmic "signal" (dashed line), and 
random pixel noise (dotted line). In order to obtain a regular 
field and suppress noise the combined spectrum is smoothed with 
a Gaussian filter; its width for this particular example is 3°. Al- 
though the noise contribution has not been completely removed, 
it has been considera bly reduced without affecting the signal too 
strongly. See Section 4.2 for a broader discussion of this issue. 



Efstathiou (1987), and smooth them with a Gaussian filter 
of 7° FWHM to model the DMR beam. White noise with 
a fixed rms fluctuation level of 3/iK is then added to this 
"cosmic signal"; this is in practice done on the pixels in 
real space, but for comparison we may also evaluate the 
contribution to the angular power spectrum, which is 



c. 



(3 M K) 2 



(18) 



independent of £. Finally, a Gaussian smoothing kernel of 
variance u> 2 , given by 



exp(-±<j 2 ^ + l)) 



(19) 



is applied to reduce the noise level, and to obtain a regular 
field. Note that the normalization factor C2 in Equation ( fTr[ ) 
is directly related to th e CMB quadrupole. As pointed out 
by Gorski ct al. (1994) its particular value determined from 
the COBE DMR sky maps is highly dependent on the spec- 
tral index n; therefore it has become common usage to quote 
the multipole Cio ~ (7(iK) 2 as a sufficiently spectrum inde- 
pendent normalization. 

So the various contributions to the angular power spec- 
trum C't for our example sum up to 



Ct + Ci 

.noise 

) ■ (20) 



In Figure [l], the contributions of "signal" and noise are 
shown separately, and combined to the full spectrum, all 
after 3° smoothing . 



With these relations and the spherical harmonics expansion 
of u, the parameters a and r may also be calculated directly 
from the angular power spectrum Ce, with the results 



a = J2(2l+l)C e , 



(16) 



3 ESTIMATING MINKOWSKI FUNCTIONALS 
OF PIXELIZED CMB SKY MAPS 

Throughout this section, we will illustrate the application 
of our method on a particular random field. In order to 
stick to a simple, analytically tractable model, we gener- 
ate a Gaussian random field. Its angular power spectrum Ct 
is chosen to reproduce the salient features of the DMR sky 
maps. Hence, we start from angular components given by 
the formula 



Ce, po 



= C 2 



r(i+y)r(2=s) 
r(£+£^)r(^) 



(17) 



3.1 Estimators for a pixelized sky map 

In order to estimate the Minkowski functionals from a dis- 
cre tized map we attempt to follow the prescription outlined 
by schmalzing & Buchert (1997) for cubic grids in three- 



dimensional Euclidean space. Note that their first approach 
based on Crofton's formula and counting of elementary cells 
is not viable since a strictly regular pixelization of the sphere 
does not exist. However, their second approach is based on 
averaging over invariants analogous to Equation (Jll]), and 
can easily be adapted to the sphere. If the random field is 
sampled at N pixels at locations x; on the sphere, we only 
need to estimate the values X 3 (xi) of the invariants from 
Equation (J12[) at each location. 



3.1.1 Covariant derivatives 

Using the well-known parametrization of the unit sphere 
through azimuth angle 1? and polar angle <p we can express 
the covariant derivatives at a point x = (1?, 93) in terms of 
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the partial derivatives^; 



W;i5 = M,i>, 

1 

U-u, = — 



Sill 17 

1 



sin ■& 
1 



cos •& 

o — 1 

sin v 
cos i? 



(21) 



+ 



sin 17 sin i? 

The partial derivatives in turn are best calculated from 
the spherical harmonics expansion 



(22) 



This is simply done by replacing the harmonic function Yi m 
with its appropriate partial derivative. Since the functions 
Yg m depend on tp via sine and cosine functions only, the 
derivatives with respect to can be obtained analytically. 
Partial derivatives with respect to i? are calculated via re- 
cursion formulae constructed by differentiating the recursion 
for the associated Legendre functi ons P™ , given for example 
by Abramowitz & Stegun (1970) 



3.1.2 Integrals over invariants 

We still have to account for the finite number of sample 
points. This is done by replacing the delta function with a 
bin of finite width A, 

5{u-v) « — l[_A/a,+A/a](it-i') I (23) 

where 1a is the indicator function of the set A, with 1a(x) = 
1 for x £ A, and 1a(x) = otherwise. The integrals sum- 
marized in Equation (hlh are then estimated by summation 

over all pixels^] 

1 N 

i=i 

For incomplete sky coverage we must restrict the average to 
the un mask ed pixels. This problem is adressed in detail in 
Section 3.2. 



3.2 Testing the estimators 

3.2.1 Complete sky coverage 

To begin with, let us look at the example without simulating 
the restrictions of incomplete sky coverage. Figure ^ shows 
the average Minkowski functionals of 1,000 realizations. 

Looking at the general features of all curves, it can be 
seen that the area vo of hot spots decreases monotonically 
from the value of one at low threshold, when the whole 

t Note that we use indices following a semicolon, such as u.j to 
denote covariant differentiation of u with respect to the coordi- 
nate i, as opposed to partial derivatives where we write indices 
following a comma, e.g. n 

§ We set the pixel weight factors Wi equal to 1, but this may be 
changed, if care is taken to preserve JZ^Li w i = L 



sphere belongs to the excursion set, to a value of zero at 
high thresholds which are not passed by any of the pixels. 
The boundary length vi starts from a value of zero for a com- 
pletely filled sphere. It reaches a maximum at intermediate 
thresholds, where the excursion set forms an interconnected 
pattern of patches and holes with a very long boundary. 
When the excursion set becomes emptier and emptier, the 
boundary length declines back to zero. For the random field 
shown in our example, the integrated geodesic curvature vi 
behaves largely similar to the Euler characteristic; the minor 
differences only become appreciable for fields with fewer fea- 
tures. Lastly, the Euler characteristic \ a ^ l° w thresholds has 
a value of two for a closed sphere. With increasing threshold, 
the Euler characteristic declines to negative values as holes 
open in the excursion set and give a negative contribution. 
This downward trend gradually stops as individual hot spots 
emerge, so a minimum develops, and the Euler characteristic 
attains positive values. Finally, more and more hot spots fall 
below the growing threshold, so their number and hence the 
Euler characteristic decreases again, reaching a final value 
of zero. 

A description of the individual curves can be found in 
the figure caption. 

3.2.2 Uncertainties through incomplete sky coverage 

In practice, a data set will suffer from incomplete sky cover- 
age. In order to estimate the uncertainties introduced solely 
by the galactic cut, we first contruct a single realization of 
the random field on the whole sky. The Minkowski function- 
als for this random field are calculated and roughly fit the 
analytical expectations, with fluctuations consistent with 
the areas shown in Figure ^. Then, we apply a series of 
straight galactic cuts with varying direction, but with con- 
stant width of 30°; this value reduces the number of pixels to 
exactly half the original value. Figure ^ shows a comparison 
of the true values for one field and the fluctuations intro- 
duced by the sample variance of the rotating cuts. Note that 
the smaller number of pixels does increase the uncertainties, 
but the average is not affected - the estimator remains un- 
biased. 



3.2.3 Boundary effects 

The previous subsection dealt with a random field that was 
first realized on the whole sky, then smoothed with a Gaus- 
sian filter, and cut afterwards. In order to determine whether 
the galactic cut affects the estimators derived above, we use 
the COBE DM R pixels and the cu stomized cut from the 
four-year data (Bennett et al. 1996). This time, we remove 



the pixels within the galactic cut before the smoothing kernel 
is applied. 

It turns out that by using this procedure, which is ac- 
tually the correct one for mimicking real data, the galactic 
cut severely affects the estimators, and leads to a systematic 
bias of as much as la. Figure ^ shows the unbiased results 
from all-sky maps already displayed in Figure ^, in com- 
parison the biased result obtained with the naively applied 
estimator. 

A straightforward procedure to remove these biases 
from the estimators is to further restrict the number of pix- 
els to the ones that lie "far away" from the cut. In order 
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Figure 2. Minkowski functionals for the example used throughout this Section. The areas indicate the averages and lcr— fluctuations of 
an ensemble of 1,000 realizations computed on all 6,144 pixels of the DMR sky cube. The analytical expectation values in Equation (jl^) 
with parameters determined from the theoretical power spectrum via Equation (|l(j) are almost exactly reproduced by the mean values 
(central lines). Fluctuations indicated by the shaded area are due to uncertainties in the parameter determination from a single map via 
Equation (p. a). They account for a large part of the overall fluctuations (empty area). 



to find them, we consider the indicator function of the cut 
itself, smooth it with the Gaussian filter, and consider the 
values at pixels outside the cut as their level of "contamina- 
tion". Now the sums from Equation ( |24| ) can be restricted 
to the pixels where the smoothed cut lies below a certain 
threshold. Figure |H] shows the results for an allowed level of 
1%; in practice, even as much as 5% produces sufficiently 
unbiased estimates. Note that while the mean values agree 
completely after applying the correction, the variance of the 
estimators has increased, simply because fewer data points 
result in poorer statistics. 

Apart from the galactic cut, point source contamination 
is another important source of incomplete sky coverage. Fig- 
ure shows the bias introduced by omitting 200 randomly 
scattered pixels. Obviously, the effect is much less pro- 
nounced compared to the realizations excluding the galactic 
cut; in fact the differences between the all-sky realizations 
and the restricted realizations are barely visible. Both the 
galactic cut and the random point cut affect roughly 3,000 
of the 6,144 DMR pixels with a contribution of 1% or above, 
so at first sight our findings appear inconsistent. However, 
they can be explained with the prominent geometric features 
- namely, almost straight edges - in the galactic cut. These 
are missing in a random point distribution, so the errors 
remain smaller and average out. 



4 EXAMPLES 




Figure 7. This map shows the earth's topography at DMR res- 
olution. The corresponding Minkowski functionals are displayed 
in Figure [j| 



4.1 The Earth 

In order to provide a familiar example that does not look 
like a Gaussian random field even at DMR resolution, Fig- 
ure ^ shows the Minkowski functionals of the earth's topog- 
raphy. The map (see Figure |^) was constructed by binning 

the Etopo5 datafljj into the DMR pixels. 



If The Etopo5 database gives elevations on a cylindrical grid 
of 5 arcminjlte snaring. The data files may he obtained from 
the net via |f tp : //walrus . wr . usgs . gov/pub/data/ ; see also Data 
Announcement 88-MGG-02, Digital relief of the Surface of the 
Earth by the NOAA, National Geophysical Data Center, Boulder, 
Colorado, 1988. 
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Figure 3. This plot shows the uncertainties due to a galactic cut. In order to separate them from fluctuations between different 
realizations in a ensemble (see Figure |i|), we chose a single all-sky realization (central line) of the random field, and applied 1,000 
straight cuts up to 30° latitude with randomly varying orientation of the equatorial plane. It turns out that the uncertainties caused 
by this deterioration of statistics are at least as strong as ensemble fluctuation. However, the correct mean value is still reproduced (the 
mean and the true value conincide in the central line of the jagged area), so the estimator remains unbiased. 



The curves reveal a number of characteristic features. 
All functionals experience a fairly sharp change at a depth 
between 6,000m and 5,000m, which is roughly the average 
depth of the seafloor. A peak of several 1,000m width and 
almost constant height of the boundary length Vi, and a 
corresponding minimum in the Euler characteristic x indi- 
cate the rise of the oceanic ridges. From 3,000m below sea 
level to slightly positive elevations, the boundary length re- 
mains largely constant, as the continental shelfs rise from the 
oceans; meanwhile, the Euler characteristic fluctuates with 
the disappearance of the oceanic ridges, and the opening of 
shallower, marginal parts of the oceans such as the Mediter- 
ranean, the Carribean sea or the Artie sea. Most of the land 
mass does not rise beyond 1,000m, so all Minkowski func- 
tionals gradually decline after this height; a few small peaks 
in the Euler characteristic may be - cautiously - identified 
with Antarctica, the Rocky Mountains, the Andes, and the 
Himalaya. 

4.2 How smoothing leads to noise reduction 

In order to obtain a regular field, and to reduce the level 
of the additive noise present in the data, it is necessary to 
apply a smoothing kernel to the data before calculating the 
Minkowski functionals. Usually, the choice of a particular 
width is largely arbitrary. Here we show the example in- 
troduced in Section ^ with different degrees of smoothing 
applied to illustrate the behaviour of Minkowski functionals 
in the presence of noise. 

The situation for 2° smoothing, where noise still makes 



an appreciable contribution, is shown in Figure O. The sur- 
face area vo is much less affected than the other Minkowski 
functionals; this is due to the fact that noise is incoherent 
and forms comparatively small hot and cold spots. However, 
these spots are almost as intense as the signal contribution, 
as can be seen from the almost equal width of all curves, 
and far more numerous - the Euler characteristic for the 
noise field alone reaches a maximum of the order of 200. 
Even though the extrema in the pure noise maps are spread 
out over the whole range of thresholds when added to the 
signal, and hence their number at a specific threshold de- 
creases, their contribution is still sufficiently high to make 
the signal appear completely different compared to the com- 
bination of signal and noise. 

In Figure |l(| where the results for 8° smoothing are 
displayed, noise is almost completely invisible in compari- 
son to the signal. Only about two dozen extrema of either 
kind (compare the extrema of the Euler characteristic) re- 
main, but since they have become extremely shallow, their 
contribution is not significant any more; the pure signal and 
the combination of signal and noise differ only marginally. 
Unfortunately, at a resolution of 8° the remaining signal 
does not carry too much cosmological information. 

With this example, the behaviour of Minkowski func- 
tionals under filtering at different scales has only been hinted 
at. The two filter widths of 2° and 8° are chosen to show 
two extremal possibilities, namely total dominance of noise 
and total reduction of noise. In practice, the intermediate 
value of 3° turns out to give good enhancement of signal, 
while preserving small-scale information as well. 
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-0.2-0.1 0.1 0.2 -0.2-0.1 0.1 0.2 

threshold [mK] threshold [mK] 



Figure 4. If the galactic cut is applied before smoothing the random field, as it should be done for real data, pixels in the vicinity of 
the galactic cut suffer from severe contamination. This leads to a visible bias in the estimated Minkowski functionals, particularly the 
integrated curvature V2 and the Euler characteristic \, although the area t>o and the circumference v\ are also affected. The shaded area 
shows average and fluctuations for the all-sky map already presented in Figure ^, while the empty area with thicker contour and central 
line indicates the same quantities for the maps biased through the cut. 



* i 




-0.2-0.1 0.1 0.2 -0.2-0.1 0.1 0.2 

threshold [mK] threshold [mK] 



Figure 5. The biases demonstrated in Figure |lj can be removed by also smoothing the cut, and excluding points where the smoothed 
cut still reaches a certain level. In the example shown, this threshold was set to a very restrictive 1%; obviously the mean values agree 
completely. In practice, a level as high as 5% might still produce reliable results. Note that fluctuations have increased in comparison 
with the results obtained from all— sky maps, as already demonstrated in Figure pi 
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-0.2-0.1 0.1 0.2 -0.2-0.1 0.1 0.2 

threshold [mK] threshold [mK] 



Figure 6. The same quantities as in Figure [j are compared in this plot, but with 200 random pixels excluded instead of the galactic 
cut. Contamination affects more pixels than in the case of a galactic cut — if a residue of 1% is allowed, 2,842 out of the original 6,144 
pixels remain, compared to 3,189 for a galactic cut. However, the estimated Minkowski functionals are more robust. 




-10000 10000 -10000 10000 

threshold [m] threshold [m] 

V 2 X 




-10000 10000 -10000 10000 

threshold [m] threshold [m] 



Figure 8. Minkowski functionals of the earth's topography from the map in Figure |jj with a 1° Gaussian filter applied. The resulting 
Minkowski functionals (jagged lines) look decidedly non— Gaussian; compare the analytical expectation values (smooth lines). Pronounced 
features in the Minkowski functional curves can be identified as corresponding to the main features of the earth's topography, such as 
oceanic ridges, continents and mountain ranges. See the main text for details. 
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Figure 9. Signal (shallow curves, empty area), noise (sharply peaked curves, empty area) and combination of both (shaded area) at 
2° smoothing. At this scale, the noise contribution contains far more features than the signal. Moreover, they are strong enough to 
persist when distributed over the whole range of the signal. Hence the combination of signal and noise displays a completely different 
morphology from the pure signal contribution. 




Figure 10. Signal (shallow curves, empty area), noise (sharply peaked curves, empty area) and combination of both (shaded area) of 
the same random field as in Figure 0, but at 8° smoothing. Although the pure noise map still features roughly two dozen extrema, their 
amplitude has been reduced to such extent that the pure signal is hardly affected when combined with the noise. 
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Figure 11. This panel shows the four-year data from the COBE 
DMR 53GHz channel, with a customized galactic cut and a 
smoothing filter of 3° width applied. Figure ha displays the corre- 
sponding Minkowski functionals calculated both from this map, 
and from the all— sky map with galactic signal dominating. 



of a stationary Gaussian random field. However, this result 
should be considered a illustration of the method rather than 
conclusive evidence, since our brief analysis probes a scale 
of roughly 4° (given by the squared sum of 2.6° pixel size 
and 3° Gaussian filter width). 

As stated above, considerable uncertainties are intro- 
duced through the estimates of the parameters fi, a and r 
entering the analytical expectation values for the Minkowski 
functionals of a Gaussian random field. In order to make this 
statement more quantitative, Figure [l3| summarizes the pa- 
rameters determined from the 1,000 mock realizations used 
for the shaded area in Figure hi Relative errors for the rel- 
evant parameters lie in the range of five to ten per cent, 
which is not too bad considering that little more than 3,000 
data points enter our analysis. 



5 SUMMARY AND OUTLOOK 
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Figure 13. Uncertainties in the parameter estimates. Param- 
eters were determined from each of the 1,000 simulated DMR 
maps, and binned into histograms. In order to give an idea of the 
amount of uncertainty, the range is extended to include zero in all 
plots. Obviously, the fluctuations arc still fairly large in compari- 
son to the mean value; to be precise, the relative rms fluctuation 
is roughly 8% for fj,, almost 6% for t, and more than 9% for yfa. 



4.3 Analysis of the COBE DMR four-year data 

As a last example, let us take a look at data that are both 
real and cosmologically relevant. Figure [ll] shows a map of 
the microwave sky as seen at 53GHz by the COB E satellite 
after four years of observing (Bennett et al. 1996). The data 



are restricted to 3,189 pixels receiving less than 1% from a 
smoothed galactic cut, when a Gaussian filter of 3° width 
is applied. Figure [li] displays the corresponding Minkowski 
functionals. Obviously, the analysis carried out on all 6,144 
pixels is dominated by galactic emission, while the field with 
the galactic cut applied is consistent with the assumption 



We have introduced Minkowski functionals of isotempera- 
ture contours as a novel tool to characterize the morphology 
of Cosmic Microwave Background sky maps. 

Using the framework of integral geometry in curved 
spaces, we were able to clarify the geometric interpretations 
of all Minkowski functionals in two dimensions. Writing all 
Minkowski functionals as spatial averages over invariants 
formed from covariant derivatives, lead to simple formulae 
for estimators applicable to pixelized sky maps, including 
a straightforward prescription for dealing with incomplete 
sky coverage. Finally, analytical formulae for all Minkowski 
functionals of the excursion sets of a Gaussian random field 
could be provided. 

The study of a simplified yet realistic model served to 
test the theoretically derived estimators in their application 
to simulated CMB sky maps. Among other tests, we checked 
whether the estimators remain unbiased when incomplete 
data is smoothed over the edge of a galactic cut, and found 
a prescription to deal with this problem while preserving as 
much information as possible. 

A number of examples provided further illustrations of 
the application of our method, and showed how to interpret 
the calculated Minkowski functionals in close touch with 
the analyzed data sets. The analysis of the earth's topog- 
raphy using Minkowski functionals explained how to cast 
the bridge from the behaviour of the Minkowski functional 
curves to outstanding features in the underlying random 
field. As a more serious application, we showed the succes- 
sive reduction of noise through Gaussian smoothing with 
increasing filter width - in the end, a complete removal of 
the noise effects from the Minkowski functionals is obtained. 
The final example briefly analyzed a COBE DMR map, with 
the not particularly surprising result that the field is consis- 
tent with a Gaussian random field on degree scales. 

Minkowski functionals combine the benefits of a sound 
mathematical framework and well-understood analytical 
possibilities with intuitive interpretations and easy appli- 
cability to real data. Hence they qualify as a method suited 
to study the Microwave sky at higher resolution, where the 
obstacle of poor statistics should not be an issue. While ex- 
periments to obtain high-resolution maps of large regions o f 
the sky are still und er development (Bersanelli et al. 199C, 
Bennett et al. 1995), testing the Minkowski functionals on 
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Figure 12. Minkowski functionals of the four-year COBE DMR map in the 53GHz channel. As expected, the full map of 6,144 pixels 
(dashed line) is severely affected by galactic emission; this results in long tails for all Minkowski functionals. Using a galactic cut, the 
functionals determined from the map (solid line) become consistent with a stationary Gaussian random field. The shaded area indicates 
the mean and variance of 1,000 realizations of a Gaussian random field with Harrison-Zel'dovich spectrum and pixel noise; normalizations 
were chosen to reproduce the parameters fi, a and t determined from the data via Equation (151). 



simulations is an important task for the future. In particu- 
lar, we need to assess their power to detect non-Gaussianity, 
and find the possibilities to estimate cosmological parame- 
ters from Minkowski functionals in an approach complemen- 
tary to the power spectrum analysis. 
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APPENDIX A: GEODESIC CURVATURE OF 
AN ISODENSITY CONTOUR 

Consider a scalar field u on a two-dimensional differentiable 
manifold 9JI. We wish to calculate the geodesic curvature of 
the isodensity contour passing through a point xp g 9Jt. To 
do this, we use a procedure outlined by |ter Haar Romcriy 



et al. (1991). fn a sufficiently small neighbourhood we can 
always find an explicit parametrization x(t) of the isodensity 
contour, with x(t = 0) = xo. The corresponding threshold is 
v = 
by 



it(xo), and therefore the contour is implicitly described 



It follows by covariant differentiation with respect to the 
parameter t thatjjf] 

u-iXi = (A2) 



must hold. So we can choose the tangent vectorr 



(A3) 



Actually this choice is not unique and reflects the freedom 
of parametrization; however, care must be taken to orient 
the tangent vector towards regions of lower values of u. Dif- 
ferentiating Equation ( |Al| ) a second time we obtain 



U-ijXiXj + u- k x k = 



(A4) 



whence we can now evaluate the geodesic curvature k of the 
isodensity contour via the well-known formula 



(A5) 



(x k x k ) 3 / 2 ' 
As our final result we obtain 

_ 2m ; iM ;2 M;12 — «flW;22 ~ U?2 u ill 
K = (u\ + W 2 2 )3/ 2 ' 



(A6) 



Note that this formula contains the covariant derivatives of 
it and therefore holds for any manifold 9Jt, regardless of the 
metric. 



APPENDIX B: AVERAGE MINKOWSKI 
FUNCTIONALS FOR A GAUSSIAN RANDOM 
FIELD 

isodensity contours of a Gaussian random field have been ex- 



cessively studied ever since the works of Doroshkevich (1970) 



on the g enus. Compre hensive overviews can be fou nd in the 
book by Adler (1981) or the famous BBKS paper (Bardeen 



et al. 198q). A highly instructive derivation of the average 



values for all Min kowski functi onals in arbitrary dimension 
can be found in ( [Tomita 1990| ). Nevertheless, we will out- 
line a calculation directly related to our approach to the 
numerical evaluation. 

A homogeneous Gaussian random field u with zero 
mean on a two-dimensional manifold 9JI is fully described 
by correlation function £(r). We wish to calculate the av- 
erage Minkowski functionals Vj(v) of an isodensity contour 
to the threshold v. Because of Equations (|ll]) and ( |i"2| ) it 
is sufficient to know the joint probability distribution of the 
field's value itself and the derivatives up to second order at 
some fixed point. 



According to Adler (1981) these six variables are jointly 
Gaussian distributed. Hence their probability distribution 
function can be written in concise form by arranging them 
into a vector u = (u, u ; i,u ; a, u ; 22, u ; i2); we have 



P(u) 



\/(27r) d det< 



: exp 



1 t -1 

-u a u 



(Bl) 



u(x.(t)) = V. 



(Al) 



II The overdot ' denotes differentiation with respect to the param- 
eter t, and summation over pairwise indices is understood. 
** 6ij is the totally antisymmetric second-rank tensor normal- 
ized to £12=1. 
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with the covariance matrix cr taken from flbmita (1990) 



(a -r -t 

r 

r 

-r v v/3 

-t v/3 v 



\ 











o \ 





•j/3 J 



(B2) 



Note that the parameters a — £(0), r = |£"(0)| and v 
£""(0) all depend on the correlation function £(r). 
Now we can perform the averages 

v (v) = J duP(u)G(u-u) 

vi(v) = J duP(u)<5(w-j/)^4 



v 2 {u) = / duP(u)5(i 



2w ; lU;2M;12 — W?!W ; 22 — W?2 U ;11 



(B3) 



by strai ghtforward int egration, and have recovered the re- 
sults of [Tomita (1990) in two dimensions 



V\\V ) 



*/2¥ exp r2?J 



(B4) 



-== exp - — 



As stated in the main text, the result depends on only two 
parameters, namely 



ff = *(0), 



A = 



2^(0) • 



(B5) 



This paper has been produced using the Royal Astronomical 
Society/Blackwell Science I^TjtX style file. 
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